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The results of the research performed under NASA 
Grant NGR-39 -005-105 consist of two publications in 
scientific journals, and a doctoral thesis. 

Dr„ Erling Nielsen completed his thesis entitled 
"Propogation of Low Energy Solar Cosmic Rays in the 
Interplanetary and Geomagnetic Field" during the period 
of this grant » Dr* Nielsen completed his work at the 
Bartol Research Foundation in September, 1974, and 
accepted an appointment at the Max Plank Ins ti tut fur 
Aeronomie, Lindau, Germany- 

Preprints of the publications entitled "Access 
of Solar Electrons to the Polar Regions" and "Angular 
Distributions of Solar Protons and Electrons" are 
appended; both papers have been accepted for publica- 
tion in Planetary and Space Science. Copies of the 
doctoral thesis are available from the Bartol Research 


Foundation. 
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Aceess of Solar Electrons to the Polar Regions 

E. Kiel. sen" and M,A» Pomerantz . 

Bartol Research Foundation of The Franklin Institute 
Swarthmore, Pennsylvania 19081 

Abstract — -The interaction between the geomagnetic and Interplan- 
etary magnetic fields is studied through Its effects upon the In- 
tensities of solar electrons reaching the polar caps during times 
® f strongly anisotropic electron fluxes in the magnetosheath, 
during the particle event of November 18, 1 968, electrons' ©f 
solar origin were observed outside the magnetopause with detectors 
aboard QG©-§, This Is the only case on record for which high res© 
lutlen directional flux observations are available for determining 
in detail the Electron angular distribution, and thus the electron 
density in the magnetosheath. 

Correlative studies of these satellite observations and con- 
current measurements by riometers and ionospheric forward scatter 
systems in both polar regions have revealed that the Initial stage 
of the associated Polar Cap Absorption event i$ attributable tq 
the prompt arrival of solar electrons. The electron flux precip- 
itating into the south polar region was equal to or larger than 
the mean directional flux in interpl anetary space, whereas over 
the north pole It was equal to or less than the backsca.ttered 
flu^» This evidence of a north-south asymmetry in the solar elec- 
tron flux at a time when the interplanetary magnetic field vector 
was nearly parallel with the ecliptic plane supports an open mag- 
netospheric model* The ratio of particle intensities in the 

* Present address: nax Planck I ns ti tut flir Aeronomie , Lintiau, 

Germany. 
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High Polar Latitude and Low Polar Latitude regions in the 
southern hemisphere is consistent with that determined at times 
when the interplanetary electron fluxes were isotropic. The 
analysis indicates that an anisotropic electron flux may be 
isotropized at the magnetopause before propagating into the 
polar regions. 
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1 o Introduction 

The common feature of all so-called open magnetospher ic 
models is that the geomagnetic field lines anchored at high 
polar latitudes are directly connected to the interplanetary 
magnetic field (imf) lines, thereby providing a direct path 
for low rigidity solar particles into the polar ionosphere 
(Dungey, 1 961 , 1 965; Axford et al. s 1 965; Mo rf 1-11 and Quenby, 
1971; Morflll and Scholer, 1972). The intensity of particles 
precipitating at each pole is related to the particle flux 
along the respective field lines in interplanetary space, and 
when this flux is anisotropic, i,e», when the intensities of 
particles propagating toward and away from the sun are unequal 
these models predict correspond i ngly different intensities in 
the two polar caps. 

Anisotropic solar proton fluxes have been used as tools 
in the investigation of the earth’s magnetic environment (Van 
Allen et al., 1971; Domingo and Page, 1971), However, similar 
observations during periods when the interplanetary electron 
flux is anisotropic have never been obtained. This is a conse 
quence of the relatively short duration (one to two hours) of 
appreciable anisotropy of solar electrons (Allum et aK, 1971) 
which precludes carrying out the required measurements of the 
electron fluxes over both polar caps with a polar orbiting 
satell i te . prior to the decay of the anisotropy. It is of con- 
siderable interest to make simul taneous . observations over the 
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polar caps because the magnetic field configuration in the open 
magnetospher i c models that have been proposed allow adiabatic 
access of the low. rigidity electrons, i n e 0 these particles 
are "good" tracers of magnetic field lines. Thus, at least 
in principle, such observations of access to the magnetosphere 
of anisotropic electron fluxes in the magnetosheath could pro- 
vide a definitive test of the validity of the open models. 

As emphasized by Paultkas (1974), observations of solar' 
particle fluxes in the magnetosphere indicate that the real 
magnetosphere is more complex than the models imply. In par- 
ticular, It appears that the solar particles have different 
modes of access to two regions of the polar cap, one charac- 
terized by open field lines [“High Polar Latitudes" ( H P L ) ] , 
the other by closed field lines ["Low Polar Latitudes" (LPL ) ] . 

In the light of observations of proton fluxes in the HPL 
regions during times of interplanetary proton anisotropies, one 
would expect to observe higher electron intensities in the HPL 
region over that pole which in an open magnetospheri c model is 
magnetically connected to the sun, and lower intensities over 
the other pole during times of interplanetary electron ani- 
sotropies, giving rise to a north-south (NS) asymmetry. 

The interaction between the incoming particles and the 
earth's atmosphere has the effect of enhancing the electron 
density in the ionosphere, which causes an increase in the at- 
tenuation of HF-radio waves propagating through the ionosphere. 
If the intensities of particles precipitating into the two HPL 
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regions differ, different attenuations result, all other factors 
being equal, hence observations of radio wave absorption pro- 
vide a direct measure of the intensities. However, if the spec- 
tra over the two HPL regions differ, or if one region is sunlit 
and the other is not, it is necessary to resort to calculations 
in order to. extract information about the intensities from the 
absorption observations. 

The solar flare of November 18, 1968, produced a particle 
event characterized by very large field aligned electron ani- 
sotropies, as observed with a satellite beyond the magnetopause 
(Nielsen et al . , 1 974 ). The direction of the imf was near the 
ecliptic plane and pointing away from the sun so that the higher 
electron intensity was expected over the southern HPL region. 

In the present paper, radio wave absorption measurements at 
high polar latitudes In the two hemispheres, obtained with 
riometer and ionospheric forward scatter apparatus, are com^ 
pared with absorption calculations based upon the electron en- 
ergy spectra derived from satellite observations in the magneto- 
sheath primarily to determine whether or not a NS-asymmetry in 
the solar electron flux occurred during this events 

2 „ Satellite Observations 

The earth-orbiting satellite 0G0-5 was located in the 
magnetosheath at a distance of 23=3 from the - earth at 1100 
UT on November 18, 1968. The Lawrence Livermore Laboratory 
electron and proton spectrometers (West et al . , 1969, 1973) 
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on board this satellite detected the onset of the electron 
event at 1045+05 UT , whereas 6 MeV protons first arrived be- 
tween 1 105-1 1 45 UT (Nielsen et al . , 1974 ). Thus the solar 
electrons arrived at the earth prior to these protons. 

Early in the event, even after proton onset, the elec- 
tron flux was/ much larger than the proton flux. In any event, 
precipitation of the higher energy protons that first reach 
the. earth will not give rise to appreciable electron density 
increases in the ionosphere because most of their energy is 
lost by ionization deep in the atmosphere {<50 km), where 
the lifetime of free electrons is short; Several authors 
(Van Allen et al , , 1964; Juday and Adams, 1969; Reid, i969) 
have used the empirical relation: /F = R • A, to relate the 

integral fluxes F of protons above some energy’E m ^ n to the 
riometer absorption A at a given frequency. R is a con-' 
stant, dependent only upon Potemra and Lanzerotti 

(1971), using solar proton data from the synchronous equa- 
torial ATS - 1 satellite, deduced R as a function of E . from 

min 

the 30 MHz riometer absorption observed at Byrd during the 
January 28, 1967 evert. At the time of interest during the 
November 18, 1968 event (1110 UT) a minimum value of the en- 
ergy of protons whi ch have arrived at the earth is E . = 46 MeV. 

The corresponding R-value from Potemra and Lanzerotti (1971) 
is R % 5,5. At this time the 46 MeV proton flux is <0.05 
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proton/cm -sec-str-MeV, and with a typical value of 3 of 
the exponent in a power energy spectrum of differential in- ■ 
tensities, we find that the proton induced riometer absorp- 
tion at 30 MHz i s A < 0 . 2 dB, which is the same order of mag- 
nitude as the uncertainty of riometer absorption measurements, 
and. an order of magnitude loss than the observed absorption 
at 1110 UT . 

On the other hand, the electrons which govern the elec- 
tron induced absorption during this event (^300 keV) ionize 
profusely at higher altitudes (^70 km) where relatively high 
electron densities can be maintained, owing to the decreased 
probability for a free electron to be attached to a neutral 
molecule or to recombine with a positive ion. Thus, at early 
times following onset, we expect the electron density in the 
ionosphere to be governed by the influx of solar electrons.- 

Comparison between measurements of isotropic interplane- 
tary electrons with electrons in the HPL regions have revealed 
that they were identical with respect to both intensity and the 
form of the energy spectrum (West and Vampola, 1971), This in- 
dicates that a solar electron propagating into the polar caps 
suffers neither energy gain nor loss. For the calculations 
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of HF-radio wave absorption that will be compared with the ob- 
servations, we will assume that the energy spectra J(E) of the 
electrons precipitating over the south pole are those deter- 
mined for the mean directional fluxes in the magnetosheath 
(curves A In Figure 1), while we assume that over the north 
pole the electron energy spectra are determined by the mean of 
the directional fluxes of particles which have been backscat- 
tered beyond the orbit of the earth and propagate toward' the 
sun (curves 6 in Figure 1). Since energy channels in the elec- 
tron spectrometer are narrow, the channel center energies have 
been used in determining the spectra. 

3 . Radio Wave Absorption Observations - 

Ground-based observations of ionospheric absorption of 
radio waves during the November 18, 1968 particle event were 
obtained with Bartol 6 s network of forward scatter systems, lo- 
cated at high latitudes. However, the only forward scatter 
1 i nk opera ti ng in the northern hemisphere at the time was suf- 
fering transmitter problems. Fortunately, the combination of 
forward scatter and riometer data have made ft possible to 
carry out this analysis, nevertheless. Some aspects of the 
absorption experiments that are relevant to the present in- 
vestigation are summarized in Table 1. 

Forward Scatter .. A forward scatter system comprises 
two stations separated by a distance of, typically, 1000 km. 
Identical highly directional antennas at each site are aimed 


at a point in the D-Iayer, midway between the transmitter and 
receiver, where the transmitted signal suffers scattering 
against electron density fluctuations, as a consequence of which 
a small fraction of the energy reaches the receiver. The scat- 
tering efficiency is proportional to the logarithm of the am- 
bient electron density at the scattering stratum, which during 
daylight at high latitudes is at an altitude of 75±2 km (Bailey 
et aK, 1 970), hence an increase in the electron density In this 
height interval leads to an enhancement in the received forward 
scatter signal. On the other hand, any increase in the ioniza- 
tion below the scattering stratum leads to greater absorption 
of the electromagnetic wave than that suffered during quiet 
times . 

Thus, during a typical solar particle event, the received 
power is determined by the difference between an enhancement 
effect caused by the electron density increase at an altitude 
of ^75 km and an absorption effect resulting from the increase 
in ionization below that altitude. 

The forward scatter records from Byrd-McMurdo (BM), repre- 
sentative of conditions in the Antarctic HPL-region, is repro- 
duced in Figure 2a. The absorption started at ^1050 UT, and 
continued to increase until 1110 -UT. The received signal in- 
tensity then remained nearly constant for about 5 minutes, until 
1115 UT when a. second onset, which is attributable to the arrival 
of solar protons, occurred. Thus, the absorption before 1115 UT 
was governed by solar electron precipitation. The absorption 
between 1110 UT and 1 1 1 5 UT was 26.5 dB; the uncertai nty i n es- 
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timating the -signal enhancement together with the limited ac- 
curacy in determining the signal level from the records com- 
bine to produce an uncertainty in the absorption observation 
of ±2 d B . 

The forward scatter absorption in the LPL-region is 
represented by Figure 2b, which shows the signal intensity in 
Byrd-South Pole link,- The absorption in the above mentioned 
time interval was 14.5±2 dB. During times of isotropic inter- 
planetary fluxes the intensity of electrons precipitating into 
the LPL region has been found to be a function of magnetic 
■local time (Vampola, 1971). We note for later use that the 
, above absorption measurements were made at ^0600 magnetic lo- 
cal time. 

Riometer . (a) - Antarctic measurements : The cosmic noise 

decrease produced by solar x-rays commenced at 1026 UT, reached 
a maximum at 1050. UT and ended at 1220 UT (Lincoln, 1970). 

Figure 3a shows the absorption measured at McMurdo, which is 
located in the HPL region. Since the absorption measured prior 
to the time of the solar radio frequency burst is induced solely 
by x-rays, as represented by -the dashed line in Figure 3a, the 
contribution must be subtracted to determine the absorption 
produced by charged particles only. Between 1110 UT and 1120 UT , 
the absorption induced by precipitating solar electrons was 2 dB, 
Thereafter, the absorption increased rapidly owing to the ar- 
rival of low energy solar protons. 

(b)- Arctic measurements .: Figure 3b shows the record ob- 
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tained with the 30 MHz riometer at Shepherd Bay/ Clearly, the 
absorption at this station before 1130 UT was <0,2 dB. The 
30 MHz riometer observations at Thule indicate onset at 1055 UT, 

0 .. 1 ± 0 » 1 dB absorption at 1120 UT and 0.2 dB at 1 300 UT (Ray 
Cormier pri vate communicat ion) , i.e., before proton onset 
the upper limit of the absorption at Thule and Shepherd Bay, 
both located inside the polar plateau, was 0.2 dB. 

Summary . On November 1 8, , 1 968, the ionospheric forward 
scatter absorption in the BM-path before 1115 UT and the riometer 
absorption at McMurdo before 1120 UT reflect the intensity of 
electrons precipitating into the southern polar plateau, i.e„ s . 
the HPL region in Antarctica which may be magnetically linked 
to interplanetary space. The difference between the above times 
arises from the ;f act that the forward scatter system is more 
sensitive than the riometer (especially to the high energy pro- 
tons, which arrive earliest). Thus, at 1110 UT the forward- 
scatter and riometer absorptions for Antarctica, which are rele- 
vant for an investigation of North-South asymmetry, ■ were 26.5 
dB and 2.0 dB, respectively. In the corresponding region in 
the northern polar cap the riometer absorption was <0.2 dB. . 

4/ - Theoretical Model of Radio Wave Absorption 

Electron Dens! ty Prof i 1 e » In order to determine the al- 
titude profile of electron density n g (h) created by the solar 
electrons, it is first necessary to calculate the ionization 
rate q(h) (number of electrons produced per cm /sec) as a func- 



tion of altitude h, 

Spencer's (1959) tabulation of the electron energy loss 
by col 1 isional ionization, in terms of the residual ranges of 
the. Incident electrons, provides the basis for this calculation 
(Nielsen, 1974). Both nuclear elastic scattering and electron 
slowing down due to collisions are taken into account. The ef- 
fect of range straggling arising from large discrete energy 
losses suffered in both radiative and inelastic collisions is 
considered to be negl 1 g i bl e , The CIRA (1965) atmospheric model, 
and collision frequencies listed by Ba i 1 ey et al « (1970), were 
adopted for carrying out the computations. 

The electrons and ions produced by the precipitating 
solar electrons will seek equilibrium with the medium in which 
they are produced. It is expedient to ignore the detailed 
chemical reactions that take place between the constituents 
of the atmosphere, and to treat the reactions collectively in 
terms of effective or macroscopic reactions represented by an 
effective recombination coefficient u e ff (LeLevier and Brans- 
comb, 1968), defined under steady state conditions as: 
c* eff (h) - q(h)/n|(h). 

We reiterate here that the purpose of this study is to 
establish whether or not there was a north-south asymmetry, 
with the greater flux over the south pole. In order to answer 
this question we will determi ne . two extrema: (1) the minimum 

flux that is consistent with the absorption measurements in 
Antarctica, and (2) the maximum flux that is consistent with 
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the absorption recorded in the northern polar area. If (1) 
exceeds (2), an asymmetry occurred. 

It is clear that the minimum value of yields the 

maximum va 1 ue of n e , and hence the maximum calculated absorp^ 
tion. Thus, by adjusting the flux of precipitating electrons 
so that the calculated and observed absorptions are equal, 
we determine the minimum value of the flux. Conversely, the 
upper limit on the fl ux corresponds to the upper limit of . a ^ 
Thus for calculations of absorptions in the southern and 
northern polar areas we will assume the lower and upper limits 
of a eff , respectively. 

In the sunlit hemisphere (in this case Antarctica),, 
where the solar ultraviolet radiation produces a high rate of 
electron detachment from negative ions, the .detachment rate is 
much higher than the ion-ion recombination rate, which implies 
that, at any given altitude, a ^ is a constant characterizing 
the atmosphere independent of the positive ion density n + . 
Figure 4 shows two altitude proviles of the daytime effective 
recombination coefficient. One (Potemra et al. , 1 969) was de- 
termined by comparing proton flux measurements in space with 
riometer and VLF observations of seven events. The other 
(Bailey et al. , 1970) was based upon auroral observations (al- 
titudes above 83 km) and measured proton energy spectra during 
the hard, early stages of a Polar Cap Absorption ( P C A ) event 
(lower altitudes) in conjunction with riometer observations. 

In the dark hemisphere (in this case the Arctic), the 
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i on- ion recombination resulting from an increase in the negative 
ion density cannot be neglected, and o^ ff is a function of the 
concentration of positive ions, i.e. } it depends upon the. in- 
tensity and energy spectrum of the precipitating solar flare 
particles. If n + increases above the quiet time level, as it 
does during a solar particle event, then a decreases, thereby, 
increasing the equilibrium . electron- density above that correspond- 
ing to values of a eff representative of the quiet ionosphere 
(Figure 4), Zmuda and Potemra (1972) determined the nighttime 
altitude profile of ot gff for the event of February 25 , 1 969* 
for which the production rate at altitudes above 62 km is larger 
than that determined in the present study for the November 18, 
1968, event. In light of the earlier comments, it now follows 
■that a eff during the November 1 8 , 1 968 event must have been.: l ess*: 
than, but perhaps close to the values valid for the February 25 
event for which the corresponding altitude profile . is 

shown in Figure 4. In our estimate of flux over the northern 
polar area we will use this a fiff profile for altitudes *62 km 
and the & e ff representing the quiet ionosphere below 62 km. 

Figure I indicates that between 1110 UT and 1120 UT, the 
precipitating solar electron flux was only slowly varying com- 
pared to the <100 sec average lifetime of the various species 
of positive ions (Bailey et al.» 1970), It takes a certain time 
for the ionosphere to reach equilibrium after exposure to a given 

ionizing agent, as represented by a time constant x - 1/(2 a -- n 

eff e 

Calculations of n g during the early, stage of the event (Nielsen, 



1974), show that at night x - 1 sec, whereas x a ] min during 
the day. We consider this to be justification for the approx- 
imation that between 1110 and 1120 UT the Ionosphere reached 
equilibrium almost instantaneously after the influx of solar 
electrons. Thus, in terms of the ionization rate and effective 
recombination coefficient, the equilibrium electron density is 
wel 1 * determi ned by an expression that is valid under steady 
state conditions. 

Absorption . Radio wave absorption is computed using 
the work. of Sen and wyller (1960), who calculated the complex 
refractive index of a weakly ionized gas in a steady magnetic 
field.. Their model takes into consideration the fact that the 
collision frequency of electrons ..with neutral particles is a 
function of electron velocity. 

Since the magnetic field lines in the polar caps are 
nearly vertical, the cosmic noise signal detected by a /riometer 
with an ideal vertical pencil beam antenna has propagated parallel 
with the magnetic field, 1 ,e. longitudinal propagation. On the 
other hand, in the forward scatter systems,. the zenith angle of 
the signal paths is nearly 85°, hence propagation is . essential l.y 
transverse. 

The calculated absorption equals that which would have 
been observed if the antennas had an ideal pencil beam gain 
pattern. A typical riometer has a broad-beam antenna, and this 
leads to a further increase in the observed absorption above 
that which would have been expected for a pencil beam antenna. 
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Assumlng that the ionization is homogeneous over the part of 
the ionosphere which affects the riometer signal , Hklund and 
Hargreaves (1968) found that the absorption increases by a 
factor M, Furthermore,, a riometer antenna responds equally 
well to the ordinary arid extraordinary components of the cosmic 
noise signal, and Bailey (1968) showed that this further in- 
creases the observed absorption by a factor -1.1 1. 

In contrast with the .riometer , the geometric properties 
of ionospheric forward scattering (the scattering takes place 
in an 'effective region around the path-midpoint with dimensions 
of the order of 100 km), make it unnecessary to apply any sub- 
stantial correction to the calculated absorption for a pencil 
beam antenna (Bail ey et al „ , .1 955) . The correction for the 
ordinary component is of the order of 1 percent and has been 
neglected (Bailey et al„, 1970). Consequently, no substantial 
correction of , the calculated absorption is required in this 
case, and the calculated and observed absorptions can be com- 
pared directly, 

5 . Comparison of. Observed and Predicted Absorptions 

Observations in the LPl and HPl Regions in the Southern- . 
Hemisphere , The absorption in the Byrd-South Pole path 

. ' ■ do 

prior to the onset of proton induced absorption was about half 
of that in the Byrd^McMurdo path Corrections for differ- 

ences in path length and radio wave frequency increase the 
ratio Ags/Agjyj to approximately 0,65, For a given form of the 
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energy spectrum the daytime absorption at the high frequencies 
utilized in this study is approximately proportional to the 
square root of the integrated flux. Assuming (1) that the 
spectral shapes in the LPL regioh ctnd.in the HPL region are the 
same (as is the case when the i nterpl anetary electron flux is 
isotropic 9 Vampola, 1971) and (2) that the cutoff rigidity at 
the midpoint of the Byrd -South Pole path is so low as to not 
significantly influence the absorption (the cutoff, 10 keV for 
electrons, is MO 5 V), then the integrated flux, and thus the 
flux in a given energy interval in the Byrd-South Pole path 
(l PL -region) is about (0. 65 ) 2 or roughly 40% of that in the 
Byrd-McMurdo path (HPl-reg i on ) . This is in very good agree- 
ment with measurements made during times of isotropic inter- 
planetary fluxes as represented In Figure 5, which shows that 
the flux of; 300 keV electrons in the LPL region (Vampola, 1971) 
was about 30% of that in the HPL region at 0600 h magnetic 
local time,- approximately the time of our observations. Thus 
we find that during times of anisotropic i nterpl anetary electron 
fluxes the ratio between the fluxes in the LPL- and HPL- region 
is comparable with that found when the interplanetary electron 
fluxes are isotropic. 

Ijorth- South Asymmetry -., (a) Antarctic The observed 
and calculated riometer and forward scatter absorptions, repre- 
senting conditions over Antarctica at early times following on- 
set of the particle event on November 18, 1968, are listed in 
Table 2 o 
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For-the riometer at McMurdo the results are in good 
agreement. The observed FS-absorpt i on at the midpoint of 
the BM-path at 1110 UT Is consistent with the value calculated 
with the of Potemra et aL (1 969), whereas the a ^ of 

Bailey etaL (1970) yields a value that is too low to ac- 
count for the observations. Thus, the Potemra-prof 11 e pro- 
vides results that agree with the observations if the calcula- 
tions are based upon spectra J(E) (Section 2) determined from 
the mean directional ..fl uxes observed at this time. 

On the other hand, for the Bailey a ^ profile to lead 
to a calculated F$ absorption -that -agrees -With , the- measure- 
ments without changing the calculated riometer absorption, an 
electron spectrum that is harder than J(£) is required. This 
would increase the ionization below the scattering height and 
decrease It above, In such a way as to increase the F$ . absorp- 
tion without affecting the predicted riometer absorption. How- 
ever, hardening of the spectrum would imply the operation of a 
mechanism for accelerating solar particles in the magnetosphere, 
a phenomenon that has not -been observed. We therefore conclude 
that, for the present purposes, the Potemra-prof 11 e is closer 
to reality than the Bailey-profile, and that J(E) is indeed the 
spectrum of the electrons precipitating over Antarctica at 1110 
UT on November 18, 1968. 

At any rate, among the reported a gff profiles that of 
Potemra yields a minimum value for the flux precipitating over 
the south polar area, and is thus the relevant one for our 


purposes . 
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, At 1120 UT the observed FS-absorption was larger than 
the calculated value for incoming electrons ■ This discrepancy 
is attributable to onset of proton precipitation at about 1115 
UT. 

The foregoing analysis leads to the conclusion that the 
electron spectrum at the top of- the atmosphere over the south-, 
ern polar cap was consistent with that of the mean directional 
flux observed in interplanetary space. In the preceding sub- 
section,, it was shown that the ratio of the fluxes in the-HPL 
region and in the L P L region, as deduced from the absorption 
observations, was in accord with that observed when the inter-, 
planetary electron fluxes were isotropic. These observations, 
at discrete geographical 1 ocations, of electron fl uxes over the 
high latitude areas are consistent with the interpretation that 
the spatial distribution of fluxes over the south, polar region 
is independent of whether the interplanetary electron fluxes 
are isotropic or anisotropic. This suggests that the aniso- 
tropic flux may have been isotropized, probably at the magneto-^. 

' i '- 1 

pause, before propagating to the south polar cap. 

■ (b) Arctic : Calculations corresponding to Arctic night- 

time conditions were based upon an effective recombination co- 
efficient representing the quiet ionosphere at altitudes <62 km, 
whereas for higher -altitudes the effects upon a e ^ of the inr 
crease in ionization were taken into account, . The predicted 
riometer absorption in the Arctic that would be produced. by 
the same incoming flux as that which fits the Antarctic riometer 
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observations at 1110 UT and 1120 UT on November 18, 1968, 
is 0.42 dB. However, the maximum value of., the observed 
riometer absorption is 0,2 dB,- It is therefore necessary 
to reduce the intensity assumed in the calculations by a 
factor < ° 0 f j > 4. 

Pi scussi on . • In open magnetospheri c models, only the 
backscattered flux has access. to the northern polar cap. The 
cal culated values of expected riometer absorption over the 
north pole that would be produced by the observed backscattered 
flux at 1110 UT and 1120 UT (Figure 1, curves B) are. 0.30 dB 
and 0,37 dB., respectively, Since the observed absorption in 
the Arctic was £0*2 dB, we conclude that the total flux pre- 
cipitating over the north pole was even less than the back- 
scattered flux. 


In the calculations of the ionization rate* backscatter- 
ing effects (i=e,, scattering of solar electrons out.of the 
atmosphere) were not specifically considered, Spencer (1959) 
did include backscatteri ng in his calculations of the dissipa- 
tion function for the plane perpendicular case, i^e,, when 
the ionizing electrons penetrate into the atmosphere normal .Ty, 
However, in our calculations and in the actual case electrons 
are incident at all angles between 0 a and 90°, and one would 
expect that the probability for an. el ectron to backscatter is 
greater when it subtends a small angle (ice,, is incident 
nearly parallel) with the top of the atmosphere than when it 
arrives, vertical ly» Because this dependence of backscatter i ng 
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upon the angle of incidence has not been taken Into account in 
the cal cul at ions representing either hem 1 sphere , our computa- 
tions probably yield too high a value for HF radio-wave ab- 
sorption produced by a given incident i sotropi c vfl ux . Thus, 
the possibility that the flux was larger than the mean direc- 
tional flux in interplanetary space cannot be ruled out, since 
this could, in principle* account for the absorption observed 
over the south pole, Consequently, adiabatic access . cannot be 
excluded. If access into the polar cap were adiabatic, the 
flux into the atmosphere over the Antarctic would be roughly 
equal to. the field aligned flux in interplanetary space (l.e*, 
about a factor 4 larger than J(E), since the particles would 
propagate in a.. ; manner such as to conserve their magnetic mo- 
ment (first adiabatic invariant). Because the magnetic field 
increases from interplanetary space to the polar caps by a 
factor 'vio^, only electrons with pitch angle 6 < 1° would 
reach the polar atmosphere. 

Zmuda and Potemra (1972) compared calculated and ob- 
served 30 MHz rlometer absorptions for 30 different cases of 
solar proton precipitation into the polar caps. From this we 
can evaluate the accuracy of the calculated absorption values. 
The mean value of the differences between the observed and 
calculated absorptions in percent of the calculated absorption 

^ A 0B$~ A CALC^ A CALC J for these 30 cases is about 30% and. the 
standard deviation of the differences is of the same order. 
Thus, an observed absorption which is -a factor between 1 and 
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K 6 times that cal .ciliated -for a given fl u-x'. 1 s consistent wi th 
the calculation within a confidence level of 66 %. Within this 
degree of uncertainty the smallest f1ux;J (E) (with a spectrum 
of the same form as J (E) that is consistent with the observed 
2 dB absorption at the south pole i s smaller than J(E), by a* 
factor (1=6)^ = 2 0 56, --i = e=, J ^ ( E ) 0 = 39 J QO On the other hand 

the largest flux J n (E) that is consistent with the upper limit 
on the absorption at the north pole is Smaller than 0{E) by a 
factor (0.42/0.20) 2 = -4.4, i.e., J N (E) « 0.23 J Q „ Thus, taking 
into account the uncertainties in the calculations, we still 
find that the flux over the southern polar cap is larger than- 
that over the northern by at least a factor 1 = 7 , 

All of these considerations lead to the conclusion that 
a north-south electron asymmetry has actually been observed for 
the first time. 


60 Con cl us ions 

Our analysis of satellite measurements of solar electron 
fluxes in the . magnetosheath during the early phase of the Novem- 
ber 18, 1968 event, together with indirect observations of the 
intensity of solar electrons precipitating into the earth's 
polar caps deduced from concurrent ground-based observations, 
has led to the following conclusions: 

1= The electron onset in i nterpl anetary space and 
the onset of High-Frequency radio wave absorption 
were essentially simultaneous, hence, the initial stage, 



- 22 ~ 


of the PCA is attributable to the arrival of solar 
electrons* 

2. A north-south asymmetry of electrons occurred* 
and the flux over the southern polar area was >1*7 
times that in the Arctic* 

3* The flux precipitating over that pole which in an 
open magnetospher i c model is magnetically connected 
to the sun (the south pole in this case) was ..equal 
to or larger than the mean unidirectional intensity 

of e] ectrons. -observed -.In -interplanetary .space. 

4.. The flux required to produce the observed ab- 
sorption over the north pole was equal to or less 
than the backscattered flux* 

5. The ratio between fluxes in the LPL- and HPL- 
region in Antarctica during a period of solar elec- 
tron anisotropy in interplanetary space is consis- 
tent with that observed when isotropy prevails. 

All of these results are consistent with the concept of 
an open magnetospheric model » and with the conclusion that an 
interplanetary anisotropic electron flux may be rendered iso- 
tropic at the magnetopause. 
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Ground-based HF-radio-wave experiments. 



Transmitter 
Locati on 

Receiver Location 

Polar Cap Region 

Symbol 

Frequency 

MHz 

Forward 

Byrd 

Me Mur do 

HPL 

BM 

23. ; 28 

Scatter 

Byrd 

South Pole 

LPL 

BS 

23,40 



Me Mur do 

HPL 


30 

Ri ometer 


Shepherd Bay 

HPL 


30 



Thule 

HPL 


30 
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Table- 2 

Observed and calculated rtome ter a nd forward scatter 'ab- 
sorption (In dB) in Antarctica on “November-1 8, 1 968 0 




FIGURE CAPTIONS 

Electron energy spectra at. early times during the 
November 18, 1968, solar cosmic ray event. Curves 
A are the spectra of the mean directional fluxes as 
observed by OGO-5 in the magnetosheath, while curves 
B are the mean spectra of electrons propagating 
toward the sun,- h e „ of the mean fl ux of -.back- 
scattered electrons. 

The intensity (in dB) vs time record for two forward 
scatter systems. When the transmitter is turned off 
for about 5 minutes every half hour* the receiver 
functions as a Home ter, recording cosmic noise, 

A Sudden Ionospheric Disturbance (SID) commenced at 
1026 UT (Lincoln, 1970), However, the average power 
received in both forward scatter systems, was constant 
between 1026 UT and the time of onset of absorption 
owing to particle precipitation, with the exception 
of- : the spike in the BS-record at 1046 UT. This in- 
dicates that the enhancement and absorption effect 
just balanced each other, and we therefore -assume 
that -the x-ray flux did not significantly affect the 
forward scatter signal during the particle event. A 
solar rf intensity enhancement lasting about 5 minutes 
commenced at 1046 UT , and this. is probably the cause 
of the peak in the . BS-signal , Only the receiving 
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antenna in the BS-link had a gain which was differ- 
ent from zero in the direction of the sun at this 
time, and that is the reason for the absence of a 
peak i n the BM-st gnal „ 

(a) Represents the pre-event reference level of the 
signal , (b) the estimated enhanced level of the sig- 
nal between 1110 UT and 1115 UT that would have been 
observed in the absence of absorption below the scat- 
tering stratum, (c) the signal Intensity in the above 
mentioned titoe interval,. The differences between 

(b) and (c) are about 26.5 -dB and 14.5 dB for the 
BM-path and BS-path, respectively, and represent 
the estimated absorption of the forward scatter-sig- 
nals. 

Fig. 3 (a) The intensity vs time records for an Antarctic 

riometer and (b) the absorption vs time profile for 
an Arctic riometer, during the event of November 
18, 1968. The three horizontal bars shown in (a) 
indicate times during which the magnetic field In 
the magnetosheath had a northward component. The 
two earliest intervals are characterized by decreases 
in the absorption, indicating that a northward field 
prohibits proton access to the magnetosphere. The 
third occurs at the time of maximum absorption, 
which might have been larger if the field had been 
directed southward . 
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Fig. 4 Altitude profiles of the effective recombination co- 
efficient a eff for daytime and nighttime conditions. 

■ Because the production rate above ^60 km during the 
event of February 25, 1969, was less than that at 
1110 UT on November 18, .1968, the -a .^-values at 
altitudes above ^60 km are larger than those expected 
on November 18, 1968. Thus utilizing a eff -values 
corresponding to the event of February 25, 1969 for 
the present purposes will lead, to too large a value 
for the calculated Intensity of electrons precipitat- 
ing over the Arctic region. 

Fig. 5 During times of isotropic solar electron fluxes in 
interplanetary space, Vampola (1971) found that the 
solar electron fluxes formed a polar plateau ( H PL 
region) with constant intensity, and at lower lati- 
tudes a quasi-trapping region (LPl region) where the 
intensity was-Jess than ,50% of that in the HPL re- 
gion and varied with magnetic local time. In the. 
figure J A is the flux perpendicular to the field 
lines in the LPL region [ J x (LPL ).] and in the HPL re- 
gion [ J A { H P L ) ] , and the curve is a 1 east-squares 
linear fit to data in which the points are the ratio 
J x ( LPL-)/ J j. (HPL ) for 31 2±27 keV electrons, plotted 
as a function of magnetic local time (Vampola, 1971). 
The data point denoted by a cross is determined in 
this study. 
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Abstract - High angular-resolution measurements of di- 
rectional fluxes o*f solar particles in space have been obtained 
with detectors aboard 0G0-5 during the cosmic ray event of No- 
vember 18, 1968. This is the only case on record for which 
sharply-defined directional observations of protons and electrons 
covering a wide rigidity range (0.3 MV to 1.5 GV) are available. 

The satellite experiment provided data for determining 

pitch angle distributions with respect to the direction of the 

local interplanetary magnetic field lines during the lengthy 

highly anisotropic phase of the event. It was found that the 

unidirectional differential intensities j (6 ) of 3- to 25-MeV 

protons varied in accordance with the relationship j ( 6 ) = 

2 1 

bQ + b-j cos e + b 2 cos 9, where bg and b-j > 0, and b 2 is posi- 
tive, zero or negative. Soon after onset, 79-266-keV electrons 
arriving from the direction of the sun displayed an anisotropic 
component with the intensity varying as cos 0. Later, afjdouble- 
peaked distribution appeared at the lower energies, whereas the 
flux at the upper end of the range covered by the experiment 
became isotropic. These results have been Interpreted in the 
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light of the temporal flux profiles and the state of the inter- 
pl ane tary med i um. 

The observation of the unusually large and long-lasting 
anisotropies lead to several conclusions including: 

(1) if injection of the solar particles was instantaneous 
the diffusion coefficient was either constant or increasing with 
distance from the sun ; 

(2) if the solar source emitted particles over an ex- 
tended period, and' there is evidence to that effect, there was 
weak scattering in the region between the sun and the earth and 
a strong scattering region beyond the earth’s orbit; 

1 ‘ • (3) solar electrons were stored near the sun; 

(4) the observed angular distribution of 200-MV protons 
in the magnetosheath was in good agreement with that dedWced 
in an earlier analysis of polar orbiting satellite observations 
and trajectory calculations. 
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1 . In troducti on / 

Satellite observations have revealed several general 
features of solar flare particle fluxes (McCracken et al., 
1970). One typical characteristic Is that the flux is gen- ; 
erally anisetrepie» In view Of the r*61§V6fiS§ §f thi iflftjliP 
distributions of these fluxes for the understanding of the 
role of the inter-pl ane-tary magnetic field in controlling the 
propagation of charged particles, an investigation of the 
pitch angle distributions in interplanetary space of low- 
energy protons and electrons of solar origin would be of great 
interest. However, owing to the lack of suitable observations, 
this has not heretofore been feasible. 

If no scattering occurred between the sun and the earth, 
the solar cosmic rays would propagate in a manner such as to 
conserve their magnetic moment (first adiabatic invariant), 
thereby producing strongly collimated fluxes at 1 AU. Since 
such strongly-col 1 imated particle fluxes have not been observed, 
their angular distribution can be considered the signature of 
the combined effects of the diverging interplanetary magnetic 
field (imf) lines and scattering mechanisms in Interplanetary 
space. Spatial fluctuations in the imf provide the governing 
scattering process (Kaiser, 1973), and thus the angular distri- 
bution yields information about the interplanetary magnetic 
field regime. S 

Theories of cosmic ray transport have been treated ex- 
tensively in the literature (Jokipil, 1968, 1971, 1972; Roelof, 
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1968, 1969 , 1974; Klimas and Sandri , 1971; Kaiser, 1973 ; Earl, 
1973a and b). We shall not discuss these in detail, but note 
that because particle transport may depend strongly on magnetic 
rigidity, it is of obvious interest to analyze observations which 
cover a wide range of rigidities. Therefore, the solar flare on 
November 18, 1968, which produced relativistic solar cosmic rays 
as well as highly-anisotropic fluxes of low-energy protons and 
electrons is of special importance, in view of the fact that 
this is the only recorded event for which well defined direc- 
tional intensity observations have been obtained for protons 
and electrons covering a wide rigidity range. 

Since this particular chromospheric eruption attracted 
widespread interest, the World Data Center has issued a special 
data compilation (WDC-A, 1970), For the present purposes, it 
suffices to note that the Importance IB flare near the west 
limb (N20°, W90°) commenced at 1026 UT, The onset of Ha emis- 
sion coincided with a rapid increase in x-ray flux level, indi- 
cating that particle acceleration apparently occurred during the 
flash phase. The i nterpl ane tary disturbance which gave rise to a 
Sudden Impulse (SI) at 1630 UT on November 18 is probably a shock 
and is relevant in the analysis of the particle fluxes. 

Although the angular distribution has been determined 
experimentally for relativistic solar cosmic rays (Figure 1) 
observed during a number of ground-level events or GLE (Duggal , 
et a!,, 1971; Maurer et a!,, 1973), no similar detailed invest- 
igations have been previously conducted for the lower-energy 
particles that are observed with spacecraft. 
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During the November 18, 1968, event, Orbiting Geophysi- 
cal Observatory 5 (OGO-5) was operating in interplanetary space, 
measuring with high angular-resolution di recti onal -fl uxes of 
low-energy protons (0.57 < T < 46 MeV) and electrons (79 < T < 
1530 keV) . 

The purpose of the present study is primarily to analyze 
OGO-5 data recorded during this event in order to determine the. 
angular distributions around the direction of the imf lines at 
early times (the first 12 hours following onset of the particle 
event), when the satellite was outside the magnetosphere. In- 
formation to aid in the interpretation of the observed angular 
distributions is obtained by studying the temporal f 1 ux |prof 1 1 e . 

2. Experimental Considerations 
A. Instrumentation 

The low-energy proton and electron data were obtained 
by the Lawrence Livermore Laboratories (LLL ) energetic-particle 
experiment on OGO-5 (West et al ... 1969, 1973). Some of the 
pertinent aspects of. the experiment that are relevant to the 
present investigation ar v e summarized here. A schematic cross 
section of the instrument is shown in Figure 2, and its relevant 
characteristics are listed in Table 1. 

In the interest of background rejection for the electron 
spectrometer, both the detectors for each energy channel (except 
the two lowest where the detectors must be fairly thick (^0.2 mm) 
to reduce electronic noise), and their corresponding background 
detectors, were designed with a thickness corresponding approxi- 
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mately to the range of an electron of the specified energy. 

Pulses representing the peak of the distribution are selected 
by pulse height analysis. This design reduces the number of 
background counts which are produced by bremss trahl ung , pene- 
trating energetic protons and degraded electrons scattered 
off the spectrometer walls. 

The proton telescope is located in the line of sight 
of the entrance aperture of the high field electron spectrometer. 

It comprises four detectors and two absorbers. The lat- 
ter are Included to Increase the energy range covered by the 
instrument. The energy of a proton is determined by measuring 
the energy deposited in each detector that It penetrates, and 
through the logical analysis of the pulses from the detectors. 

The proton channels and the logic statement discriminator set- 

•• . Y 

tings used in defining the energy channels pertinent to this 
investigation are specified in Table 1. Note in the discussion 
that proton channels P3, P4, P5 and P6 are often referred to in 
terms of their approximate mean energy 1, 3, 9 and 25 MeV. The 
instrument had two additional channels extending down to 100 
keV. However, protons of non-solar origin in this energy range 
are often found in the magnetosheath and nearby upstream-wave 

o 

region (.Lin et al., 1 974; West and Buck,. 1974) and hence were 
eliminated from the present investigation. 

The intense low-energy portion of the electron spectrum 
(<4 MeV) is prevented from reaching the telescope by the magnetic 
field in the electron spectrometer; also 9 ^the detectors are so 
thin that, usually, an electron cannot deposit enough energy 
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in the detectors to overcome the thresholds. The same argu- 
ment holds for bremss trah 1 ung . The primary source of back- 
ground is penetrating energetic protons. Background correc- 
tions for all of the particle detectors are made according to 
the procedure outlined by West (1972). 

OGO-5 was sun-earth oriented, hence the experiment aper- 
ture had to be scanned relative to the stabilized spacecraft 
to obtain directional information. NASA provided a special 
scan platform mounted on one of the spacecraft booms which always 
pointed toward the’ earth; hence the experiment scanned so as to 
look perpendi cu 1 ar to the earth-satellite radius vector (note 
the bottom panel of Figure 3). The scan was 115 degrees at 
a normal rate of 3° per second but occasionally 1 1/2° per second 
(note that the longest data acquisition interval was 4.6 seconds 
so that the angular resolution was not too adversely affected 
by the scan rates). It is unfortunate, certainly for this in- 

i ’ / 4 

vestigation, that a more complete (+180 degrees) scan was not 
avai 1 abl e . 

Simultaneous magnetic field data were obtained with 
the UCLA triaxial fluxgate magnetometer (P. J. Coleman and C. T. 
Russell) on board OGO-5. Because the primary purposes of the 
data analyses reported here are to determine the angular dis- 
tributions and to deduce their physical significance, it is 
essential that the magnetic field and particle observations 
were conducted simultaneously. 4^ 

Maximum information about pitch-angle distributions 
can be deduced from measurements recorded during periods in 
which the plane of scan contained or nearly contained the direc- 
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tion of the magnetic field lines. On the other hand, when the 
scan was perpendicular to the field lines, only observations 
near 90° pitch angle were obtained. To exploit the instrument 
to the fullest extent for the present purposes, the satellite 
must be outside the magnetosphere at'a position where the plane 
of scan contains the garden-hose field line. During the 3 1/2 
year lifetime of 0G0-5, these conditions were satisfied only dur- 
ing the GLE on November 18, 1968, 

The relevant part of the orbit of 0G0-5 from 0200 UT 
on November 18, 1968 to 1500 UT on November 19, 1968, is shown 
in the upper panel of Figure 3,' The bottom panel shows how the 
experiment aperture was scanning at the start of the measure- 
ments, Figure 4 shows the time profile of the magnitude, and 
longitude and latitude in Geocentric Solar Ecliptic (GSE) coordi- 
nates, of the observed magnetic field vector. These observations 
are used to determine the location of the bow shock and the mag- 
netopause, 

B, State of the Interplanetary Medium 

During this event, Durney et al. (1972 ) and Qiienby et al, 
(1974) observed the interplanetary magnetic field on the sunward 
side of, the bow shock with Explorer 35. Their measurements indi- 
cated a rapid change in the direction of the imf at about 1630 UT, 
following which there was a decrease in the magnetic field in- 
tensity at 0G0-5 which was in the magnetosheath during the inter- 
val 1600-2100 UT (Figure 4), At the same time (1630 UT) the 
geomagnetic field was perturbed, probably by a magnetohydrody- 
namic shock propagating in the solar wind, giving rise to a 
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Sudden Impulse (SI) recorded at low-latitude magnetic observa- 
tories (Kawasaki and Akasofu, 1970) 0 After about 1800 UT the 
imf -remained nearly aligned with the sun-earth line until 2400 UT. 

■ We have indicated in Figure 3 the direction of the field 
in the magnetosheath which is tangential to the magneto- 
pause. All points to the east of this line are magnetically 
connected to the sun, whereas all between it and the dusk and 
evening magnetopause are not. Also, particles propagating from 
the sun with trajectories on the sunward side of this tangen- 
tial field line do not interact with the magnetosphere. 

On the sunward side of the magnetic irregularity, be- 
tween 2200-2400 UT, the average longitude of -.IT is 55°, as in- 
dicated by in Figure 3. 

To illustrate the directions in which fluxes are meas- 
ured, the look directions of the particle detectors are shown 
in Figure 5. 

3 . Intensity Variations 

A. Measurements 

Before presenting the observed temporal variations of the 
particle intensities and angular distributions a few points should 
be noted. It can be shown that adiabatic deceleration causes 
a maximum change in rigidity in the time interval of interest 
( r v/|2 hours) which is less than 15% (Nielsen, 1 974); hence, 
the effect of deceleration is insignificant and is not a cause 
of the phenomena to be discussed. Owing to the form of the 
energy spectra and the low velocity of the solar corpuscular 
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stream in the magnetosheath, it is evident that the corrections 
for convective effects (Forman, 1970; Balogh et al . , 1973) are 
insignificant because of the large diffusive anisotropies en- 
countered in this study; hence they have been ignored. 

Figures 6 and 7 show J, the observed particle flux aver- 
aged over the plane of scan, plotted as a function of time for 
electrons and protons, respecti vely . In all figures displaying 
data points, the error bars represent twice the standard devia- 
tion of the mean, determined from a variance analysis, taking 
into account time variations of the flux. The intensity of 9- 

and 25-MeV protons and of the electrons increases rapidly to 
a maximum and then starts to decrease. At the time-of arrival 
of the magneti c~1rregul arity there is a decrease in counting 

rate in all the electron channels and in the 3- and 9-MeM ; ^ 

proton channels, while the 25-MeV channel is affected less. 

There is a small maximum in the 79-keV electron flux at 1530 
UT, while no indication of a maximum is observed at higher 
energies . 

B . Analysi s 

Coherent Propagation . Earl (1973b) found that wiien the 
imf power spectrum is steep, a pulse, of particles from a solar 
flare (the coherent pulse) may arrive at the orbit of the earth 
preceding the diffusing parti cl es s thereby creatingan initial flux 
maximum in the intensity profile. The Earl theory also treats the 
arrival of the initial group of particles predicting a propagation 
speed along the spiral interplanetary field >80% of the free propa- 
gation speed. Figure 8 shows the time-dependent arrival of the 
first particle fluxes in the various energy channels. The. straight 
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line with slope-1 which best fits the data yields a total path 
of 1.74 AU. The Earl theory indicates a field-line length 
of M . 4 AU which is in accord with previous observations; Lin 
(1974 ) has reported M.2 to 1.5 AU . . 

The maximum in the distribution (see Figures 6 and 7) 
predicted by coherent propagation shou 1 d appear at t^ £ 1.25/3$ hr. 
following particle injection (here we have assumed a distance 
traveled of 1.5 AU , an upper 1 imi t of expectations ) . The observed 
time of maximum of 79-keV electrons is not inconsistent with 


t+ie Earl theory, and hence the initial peak in this, case would 
appear to be the result of scatter-free coherent propagation. 

For the other cases, excluding P3 for which the analysij|is con- 
clusive, the calculated times of maximum are too early to 
support a picture of coherent propagation alone. 

Diffusive Propagation . Applying the theory of anisotropic; 
diffusion with a boundary (ADB) of Burlaga (1967) to the proton 
measurements, as shown in Figure 9, we find that the data representing 
the intensity increase show a good fit to a straight line; this is 
in. accordance with the theoretical prediction. ' We further find that 
the data corresponding to the early stages of the decay, similarly, 
are in accord with this»model. From the slope of the curves in 
Figure 9, the time of maximum is determined to be 1513 UT and 
1327 UT for 9- and 25-MeV protons, respectively, consistent with 
observations (Figure 7), and 1718 UT for 3-MeV protons. The 
smooth intensity profiles corresponding to the straight lines 
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in Figure 9 are plotted in Figure 7, to highlight the proton 
intensity variations caused by the magnetic irregularity (see 

Figure 4). Owing to short rise times and small Intensity 1n- 

\ 

creases, a similar analysis of the electron and neutron monitor* 
data was inconclusive. . ' 

The mean free paths, calculated using the Green func- 
tion solution to the diffusion equation for a spherically sym- 
metric system with a point source at the sun, are found to be 
X ii = 0.12, 0.11 and 0.12 AU for 3-, 9- and 25-MeV protons, re- 
spectively. 

Associated Magnetic Field Intensity Variations . The 
sudden decrease in the proton and electron fluxes at 1530 UT 
(Figures 6 and 7) coincide with the sudden decrease in the mag- 
netic field Intensity, and Is most likely caused by It. Because 
the magnetic flux through a cross-sectional area of a magnetic 
flux tube is constant (V*Bv=0), the observed decrease in the 
field intensity by a factor 2 (Figure 4) must produce a corres- 
ponding increase of the area, leading to a reduction of the par- 
ticle flux inside the magnetic irregularity to half the flux out- 
side. This is in accor'd with the observations. The decrease in 
J inside the magnetic irregularity is to some extent caused by 
the fact that only particles with pitch angles near 90° are de- 
tected there. Thus, owing to the restricted angular coverage, 
the angular distribution is undetermined. 

Summary . ~~It has been shown that the first peak in the 
proton flu* can be described in terms of diffusion theo.]|y . If 



the protons did diffuse from the sun toward the earth, the 
diffusion approximation, i.e» a first order approximation 
valid for small anisotropies, probably yields a realistic 
value of the mean free path which we have determined to be 
^ ^ 0.1 AU. Early in an event the anisotropy is predicted 

to be large and to decay inversely with time (Fisk and Axford, 
1969) . In the early phase of the event under investigation 
the anisotropy remained large, ^100% (Section 4), and this may 
raise doubts about the validity of results obtained using the 
diffusion approximation. Early in an event the contribution 
of backscattered particles to the density is small, and thus 
of no governing influence upon the observed density. Absence 
of backscattered particles-^i . e . , persistent large anisotropies-- 
possibly owing to a diffusion coefficient rapidly increasing 
with heliocentric distance, would therefore not seriously affect 
the cosmic ray density. If the cosmic rays in fact are diffus- 
ing between the sun and the earth the observed density would 
thus fit the diffusion equation to yield a mean free path 
which is a good approximation of the conditions at helio- 
centric distances < 1 AU. 

Furthermore, -the time of maximum flux in the 79-keV 
electron channel and the time of arrival of these electrons 
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are in line with expectations based upon the theory of co- 
herent propagation . 

It was also found that the intensity of 79-keV elec- 
trons exh i bi ts a smal 1 secondary maximum at about 1530 UT, 
i.e., coincident with the detection of. the magnetic irreg- 
ularity,. We will discuss these results later after pre- 
senting the angular distributions in the next section, 

4, Angular Distributions 


A . Measurements 

The angular flux distributions around the direction of 
the observed magnetic field of both low-energy protons and elec- 
trons have been determined as they change owing to spatial and 
temporal effects. 

The shock propagating from the solar direction stretches 
out the magnetic field lines in the region of space it has tra- 
versed causing weaker magnetic scattering behind compared with 
the region ahead, and this may account in part for the observed 
long duration of the large proton and electron an i sotropies * 



Protons . Representative angular distributions of pro- 
tons are shown in Figure 10. At 1335 UT (i.e. at a distance of 
700 Earth Radii (R E ) ahead of the magnetic discontinuity) the 
intensities of 3- to 25-MeV protons vary over the angular inter- 
val where the flux is measurable according to 

j(p) « b o +b^y+b 2 u 2 . (1) 

A statistical analysis indicates, with a 95S» confidence level, 
that for the 3-MeV proton flux h 2 =0, while at 25 MeV b 2 <0. 
Thus, the angular distribution is energy dependent in the inter 
val 3- to 25-MeV. Owing to the lack of information about the 
magnetic field configuration in the magnetosheath, a theoretical 
Investigation of the distributions is not feasible, 

Durney et al . (1972) utilized observations of 100- MeV 
. (400 MV) protons made by'.a. polar orbiting satell ite together 
with trajectory calculations to deduce their angular distribu- - 
tion in interplanetary space at early times during the event 
of November 18, 1968. This result is shown in Figure 11 tor 
gether with our observation of the distribution of 25- MeV 
(200. MV) protons. The angular distributions of protons with 
rigidities of 200 MV and 400 MV, both represent conditions 
shortly after. onset of. the event, and in view of the small 
difference in rigidities, the good agreement indicated in the 
figure constitutes heretofore unavailable evidence supporting 
the validity of the trajectory calculations. 

Figure 12 shows the proton angular distributions ob- 
served within ^300 R e earthward of the magnetic discontinuity. 
The most striking features are the minima near 90* * with a 
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linear distribution for positive y and a small maximum in 
the distribution for y - -0.4. This pattern is clearly dis- 
played by 25-MeV protons, but is barely discernible in the 
♦ 

distributions of 3- to 9-MeV protons. The angular coverage is 
limited inside the magnetic irregularity, hence the results 
there are inconclusive. 

An angular distribution representative of the observa- 
tions on the sunward side of the tangential discontinuity for 
all energies is .shown in Figure 13a. It is consistent with an 
exponential representation. However, owing to the statistical 
spread of the data points only a second order term in an ex- 
pansion in cos 6 can be determined. ,j.;< 

At *2230 UT for 25 MeV and at *2230 UT for 3- and 9-MeV 
protons the angular distributions change Into a superposition of 
an isotropic component upon an anisotropic component which is 
linear in u for y > 0 only (Figure 13b). At the same time the 
intensity of the isotropic component starts to increase rapidly, 
leading to a reduction of the anisotropy as the satellite 
approaches the magnetopause, and culminating in a nearly iso-' 
tropic proton distribution in the magnetosphere near the mag^ 
ne topause . 

El ectrons . Representative distributions of 79-keV elec* 
trons are shown in Figure 14. where (a ) is the angular variation 
in directional fluxes at early times, i.e. <45 minutes following 

onset, and (b) illustrates the features at later times, In 

- , 

both cases the Intensities for ji > 0 are linear In cos 9. A 
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statist i cal analysis reveals that, while the distribution of 
the backsca ttered flux (p < 0) at later times (b) is independent 
of 0, at early times (a) the intensity decreases with cos 9. 
•Similar results are obtained for the 158-keV and 266-keV 
electron flux, although in these cases the transition to a 
constant backscattered flux occurs earlier. 

The temporal variation of the electron anisotropy is 
shown in Figure 15. For 79-keV electrons it is characterized 
by a gradual decrease with time, while at higher energies there 
is a pattern of an initial fast decrease (during 100 minutes 
after onset) followed by a nearly constant anisotropy until 
the arrival of the magnetic irregularity, when the anisotropy 
becomes indeterminate owing to the limited angular coverage. 

Figure 16 shows several representative distributions 
jj 1 •/ obtained during the time interval 2200-2400 UT, when (1) the 
satellite is > 3 R £ from the magnetopause, (2) magnetic field 
measurements are available and (3) the angular coverage is 
adequate. Figure 16a reveals a statistically significant bi- 
directional anisotropy in the 79-keV electron flux. The 158- 
keV electrons (b) seem 1 to exhibit the same behavior though 
less pronounced, while at 266 keV the electron flux is iso- 
tropic. 

B, Di scussi on ; 

Protons . The observed proton distributions inescapably 
exhibit one common feature, namely that the anisotropy invariably 
is very large. The solar protons appear to be steadily streaming 

past the orbit of the earth; but even after this has continued 

! ' . 
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for 12 hours, hardly any flux propagating in the opposite direc- 
tion .along the field is detected. 

The analysis .shows that at least 4 hours elapse after 

onset before the proton flux from the antisolar direction rises 

• , ft. ' 

above background. Hence the diffusion mean free path, 'X„ , of 
protons propagating parallel to the direction of the mean imf 
beyond the orbit of the earth, appears to be at least half the 
distance traveled in 4 hours, i.e. X„ for 3-, 9- and 25-MeV de- 
rived this way is ^1.2, 1 .9 and 3.4 AU, respectively. McCracken 
et al. (1967) using similar arguments, found for some events 
X„-values up to 2.7 AU for 10-MeV protons. However, Jok-1.pt i- . 
(1968) argued that an estimate of X« on this basis leads to 
an erroneous result. 

In the theory it Is envisaged that any net change in 
pitch angle is a sum of many small angle changes, and the time 
required to produce a backscattered flux propagating upstream 
along the field lines thus depends critically on the variation 
of scattering rate with pitch angle. Because <(Ap) >/At is 
small around 6*90° it may take an exceedingly long time for 

i ■ 

particles to be scattered through 90°, thereby giving rise to 
large backs catteri ng times (n,4 hours as mentioned above). 
However, the diffusion mean free path X»« is governed by the 
scattering in the region 0 < e < 90°, and thus it is misleading 
to directly associate Am with the backscatteri ng time. 

Jokipii finds that the mean free path, cal cul ated on the 
basis of the backscatteri ng time, is roughly an order of magni- 
tude larger than the diffusion mean free path deduced from the 
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intensity-vs-time profile. Thus, we f i rid \ ■v 0(0.1 AU). 

II 

Jones et aK (1973) found that the diffusion coefficient 
near 0=90° is substantially different from zero, giving rise 
to a considerable sc attering rate. In view of the controversy 
over the rate of scattering of particles with pitch angle near 
90 , we conclude that the magnitude of the diffusion mean free 
path of protons beyond the orbit of the earth lies between the 
approximate limits 0„1 AU and 1 AU„ The mean free path during 
this. event was calculated for 1 MeV protons by Quenby et al » 
(1974) using magnetic field power spectra. They find X deter-. 

it 

mined in this manner to be a factor ^12 smaller than the value 
deduced from the observed time-lapse between event onset and 
rise above background of the flux from the antisolar direc- . 
tion, a result that is in line with our findings. 

— e .ctron s . The solar-electron angular distributions 
exhibit two dominant features. Superposed upon an isotropic 
flux is an anisotropic component with pitch angles that do 
not exceed 90°, and for which the intensity is proportional 
to cos 0 . Thus, while part of the electron population isiso- 
tropic, there is also a steady flow of electrons past the orbit 
of the earth, a situation which is not characteristic of dif- 
fusive propagation. 
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One possible explanation of the form of this distribu- 
tion is that the electrons are weakly scattered in the region 
between the sun and the earth, whereas at larger distances from 
the sun the scattering probability increases, giving rise to 
the isotropic component However, it is not necessary to in- 
voke the effects of a distant scattering region in order to 
interpret these 'observations ; one may, instead, ascribe the 
isotropic component to large-angle scattering of the elec- 
trons out of the collimated beam. 
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Earl (1973a) found that the angular distribution in 
the case of coherent propagation should be relatively uniform 
over the forward hemisphere. However, the effect of the diverging 
nature of the interplanetary field was not included in the cal- 
culations, and thus an observed distribution might be expected 
to be more collimated than the theory predicts. It should be 
noted that while the temporal profile of only the 79-keV flux 
could be interpreted in terms of coherent propagation, the angu- 
lar distribution at that energy is similar to that observed at 
higher energies. Hence, the form of the measured angular dis- 
tribution, per se, can not be regarded as supporting the coherent 
propagation model . 

The initial rapid decrease of the electron anisotropy 6 
(Figure 15) is consistent with observations by Allum et al. 

(1971), who reported the results of a study of the anisotropies 
of electron fluxes during eighteen prompt so 1 ar-el ec tronf§bvents 
that were observed by Explorer 35. They found that the ani- 
sotropies of >_ 70-keV electrons were field aligned initially, 
with amplitudes generally ranging from 30-60% (the largest ob- 
served was 85%, July 30, 1967) that decayed to £10% within an 
hour or two. The fast decrease of 6 for 158- and 266-keV electrons 
is a typical temporal feature of prompt electron events. However, 
the initial anisotropy in the present event is unusually large 
(90 to 100%) and is regarded as evidence of weak scattering. 

Lin (1974) has noted that anisotropy measurements aboard 
earth-orbiting spacecraft may not reflect interpl anetary con- 
ditions since M0 to 10 Hz waves of terrestrial origin (Fairfield 
1 969 ; Russell et al. , 1971) are commonly observed far upstream 



from the earth's bow shocks. Since these waves seem to be assoc- 
iated with imf lines that connect to the bow shock they are far 
more probably on the dawn side of the earth than on the dusk side 
It is not likely that these waves affected the present results, 
but If they are efficient in scattering the low-energy electrons, 
we observed, this would most likely imply that the angular dis- 
tributions are more collimated in interplanetary space than at 
the point of observation. Thus, confirmation of this effect 
would lend further evidence to the concept of weak scattering 
in interplanetary space. 

Figure 15 shows that the anisotropy does not die away 
within one to two hours as Allum et al. (1971) found to^be char- 
acteristic in the usual case. However, these authors reported 
observations made during one event, on October 30 , 1967 , which 
deviated from the typical pattern in a manner similar to the 
present one. That event displayed the following features: (1) a 

double-peak structure with a minimum at almost 0330 UT appeared 
in the temporal flux profile following onset at ^0030 UT; 

(2) the anisotropy exceeded > 30 % for three hours following the 
onset of the electron increase; (3) a sudden impulse (SI) was 
observed at 0426 UT (Solar and Geophysical Data, 1967). The 
variation in^ the H-component as a typical low-latitude, magne- 
tometer station (Tangerang) was +35y for’ the November 18 event, i 
whereas on October 30, 1967, the Si-magnitude was +45y. Thus, 
these two electron events share three features in common which 
may have their origin in the presence of a hydromagnetic shock 
in interplanetary space that gives rise to the sudden impulse. 
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At later times in the 0G0-5 data, a second maximum in 
the angular distribution appears, denoting a flow of electrons 
from the antisolar direction centered on the magnetic field 
vector. The 79-keV electron flux is observed to remain ani- 
sotropic for about 12 hours. 

Bi-directional anisotropies have previously been ob- 
served only for low-energy protons in so-called Energetic Storm 
Particle (ESP) events (Rao et a!., 1967), and it is therefore 
of interest to see if the present bi-directional electron ani- 

y 

sotropy can be interpreted in these terms. Such an event con- 
sists of particles swept up by a blastwave propagating in inter- 
planetary space. Some of the particles leak'out of the blastwave 
to produce anisotropic fluxes propagating both toward and away 
from the sun „ 

The gyroradius of a 260-keV electron in a 5-y field 
is ^0.08 R^. The radius of curvature of the imf lines ob- 
served around 1630 UT is approximately 20 R £ , so in the absence 
of small-scale magnetic irregularities, the electrons could pass 
adiabatically through the “kink". Though a small maximum in 
the 79-keV electron flux is observed near the time of the ar- 
rival of the shock, the profile in this case does not have 
the characteristics determined for proton ESP-events, namely 
an abruptly increasing and decreasing intensity phase and a to- 
tal duration of approximately six hours. Furthermore, we might 
expect a blastwave to be most effective in “sweeping up" par- 
ticles which are slow compared to its speed, in which case the 
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characteristics of an ESP-event in the detected low-velocity 
protons as well as in the electron flux would be observed. 
However, no proton intensity Increases occurred at this time 

V 

(around 1630 UT). This suggests that the Interpretation of 
the bi-directional anisotropy of the low-energy electrons as 
a manifestation of the mechanism discussed by Rao et al. (1967) 
is not applicable in the present case. 

The maximum from the solar direction in the observed 
bi-directional anisotropies, may indicate that the flux of 
79-keV solar electrons remained anisotropic throughout the 
entire period covered by our observati ons ( 1 ,e 'vl 4 hours fol- 
lowing onset of the flare). 

.; ; £ 

To account for these observations, which are In sharp 
contrast with the general pattern of prompt solar-electron events 
as reported by Allum et al. (1971), it appears to be necessary 
to assume storage of the observed low-energy electrons near the 
sun ( Simnett , 1 971 , 1973). Since 266-keV electrons become iso- 
tropic before 79-keV electrons, the storage time for the latter 
i s 1 onger . 

i ,• 

Summary . Considering the backscatter time we estimate 
that the mean free path of protons propagating beyond the orbit 
of earth lies between the approximate limits of 0.1 and 1 All. 

The long-lasting large anisotropies can be seen, at least in 
part, as a consequence of the weak magnetic scattering in inter- 
planetary space. Even weak scattering will eventually isotropize 
the flux, hence it is not necessary to Invoke any additional 
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mechanisms in order to account for the fact that the solar 
cosmic rays eventually become isotropic* However, we should 
note that our observations are also qualitatively eons i stent 
with the presence of a diffusive region located at a distance 
beyond the orbit of the earth (Roelof, 1974). In this model, 
cosmic rays suffer only little or no scattering between the 
sun and the earth, and diffusion does not set in until the 
particles reach the distant scattering region, where the iso- 
tropic component near the earth's orbit has its origin. The 
electron observations are in good accord with the predictions 
of this model ; the proton results, although inconclusive, 
are not inconsistent. To maintain long-lasting large field- 
aligned proton anisotropies at the orbit of the earth in the 
case of weak scattering in the inner solar system it would, 
however, be necessary for the solar-proton source to remain 
active over an extended time interval. Otherwise, in the case 
of both instantaneous injection and scatter-free ^ propagati on , 
the space within the earth's orbit would soon be emptied of 
particles propagating directly from the sun. 

i • 

5 . Cone! usions 
We have presented detailed observations of the angular 
distributions of low-energy protons and electrons during the 
•early phase of the solar-parti cle event of November 18, 1968. 
The results and conclusions are as follows: 

1) The protons displayed very large (nearly 100%) 
anisotropy, and the electron anisotropy was unusually large 
‘ arid Tong lasting. 
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2) The mean free path A„ of 3-» 9- and 25-MeV protons 
propagating beyond the orbit of the earth was between ^0,1 All . 
and 1 AU ; A„ for protons propagating from the sun to the earth, 
as deduced from the time profile of the observed intensity, was 
'V/O.l AU . For instantaneous in jeetion* this indicates t ^a t the 
diffusion coefficient was either constant or increasing wi th 
distance from the sun. 

3) If it is assumed that the solar source was emitting 

% 

particles over an extended period (i.e. not instantaneous injec- 
tion), then our observations of angular distributions are con- 
sistent with the existence of weak scattering in the region 
between the sun and the earth and bf a strong-scattering region 
located beyond the orbit of the earth. 

4) From the unusually long duration of electron ani- 
sotropies we have inferred that the solar electrons were stored 

, - o , 1 

near the sun . 

5) A field-aligned bi-directional anisotropy in the 79- 
keV electron flux was observed at late times. This phenomenon 
may be associated with ^the proximity of the magnetosphere at the 
time of observati on . 

6) Good agreement was found between the angular distribu- 
tion of 200-MV (^25 MeV) protons observed in the magnetosheath 

and that of 400-MV (MOO MeV) protons deduced from observations 
made by a polar-orbiting satellite together with trajectory cal- 
culations (Durney et al . 8 1972), supporting the validity of these 
trajectory calculations., 
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Table l: Relevant characteristics of Lawrence Livermore Laboratories energe t i c -part i c 1 e experiment 

aboard OGO-5. The logic statement is exemplified as follows: the expression 

D i (1.35) [D~i"( 5 . 4 ) +D 2 ( 0 * 5 3 ) ] requires for an output that the signal in D-j be greater than 
1.35 MeV and , i n the brackets that the signal in D ^ be not greater than 5.4 MeV o_r the 
signal in D 2 be not greater than 0.53 MeV. 


Channel 

symbol 


Energy 


Geometric 

factor 

cm 2 -keV- 

str 


Logic -statement discriminator settings 



0.57-1.35 MeV 
1.35-5.40 
5.60-13.3 
14.0-46 


1 .30x10 
1.30x10 
1 .25x10 
1.72x10 


D t (0. 57 ).ED' 1 (Y .35) +D^ 
D 1 (1 . 3 5 ) [ D 1 ( 5 . 4 ) +D^T 
D 1 (1 .35)D 2 (1 .3)[D^ 
D 2 ( 0 . 5)0 3 (1 . 2 ) [D~C5“ 


Omnidirectional 

proton 

measurement 


>80 MeV 
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Figure Captions 

Fig. 1 Relativistic solar cosmic ray intensity measured with 

ground-based nucleonic intensity detectors. Results are 
expressed as percent of the galactic cosmic-ray back- 
ground during an early epoch of the November 18, 1968 
event and are plotted as a function of the cosine of the 
angle between, the asymptotic direction of viewing (i.e., 
the direction from which the cosmic rays, causing an. in-, 
tensity increase at a given station, arrive at the boun- 
dary of the geomagnetic field) and the axis of symmetry 
in the angular flux distribution + is the centroid of 
data points, and the straight line connects it to the 
origin. The axis of symmetry is approximately paral- 
lei with the direction of the lines of force of the 
interplanetary magnetic field at the orbit of the earth h 
(after Maurer et al . , 1973). In the present paper, a 
similar analysis of low-energy-particle data, obtained 
in space with OGO-5 during the same event is carried out. 
Fig. 2 The high-energy electron spectrometer and proton tele- 
scope. Here we ^show only one of the 180° first-order- 
focusing spectrometers comprising the el ectron-detection 
. system (a second spectrometer provides four 1 ow ^energy 

channels). Note that E65 » EB6 and EB7 are background 
channels. The proton telescope takes advantage of the 
electron spectrometer magnet to completely remove electron 
background . \ 
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Fig. 3 Experimental conditions. The Geocentric Solar Ecliptic 
(GSE) coordinate system is used in which points 

toward the sun along the sun-earth line, and the Z^^-axis 
points toward ecliptic north, Y^^ i s determined by the 
right hand rule. A part of OGO-5's trajectory on Novem- 
ber 18-19, 1968, is constructed, by rotating each point 
along the orbit into the ecliptic plane (i.e. f the XY- 
plane) around the X-axis. The earth's bow shock and 
magnetopause are drawn through the points of observation 
• at 0200 UT on November 18 and 0300 UT on November 19, 

respectively: V* is the solar-wind velocity vector; 

its deflection upon crossing the bow shock is illustrated. 

B 1 and B 1 1 show the direction of the observed magnetic 

m - - m 

field averaged over the time intervals 1100-1530 UT and. 
2200-2400 UT, respectively. The bottom panel illustrates 
the satellite observational configuration; the elevation 
of which is 35°N at 1100 UT (November 18) and 21°N at 
0300 UT (November 19). 

Fig.. 4 Magnetic field observations by 0G0-5 from 1100 UT on 
November 18 to 0530 UT on November 19, 1968. 0 min is 
the minimum angle between any observed directional flux 
and the direction of the magnetic field. Between 1600 
and 2130 UT, the angular coverage extends only to about 
+10° around the direction perpendicular to the magnetic 
field vector; this is insufficient for the determination 
i of angular distributions. The flux of 3-MeV protons 

averaged over the plane of scan is also shown. 
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Fig. 5 The projections on the ecliptic plane of the directions 

in which fluxes are observed representative of the time 

interval (a) 1100-1700 UJ and (b) 1 700-2400 

The angle between the direction of observation 

and the ecliptic plane is indicated in each case. 

Fig. 6 J, the directional electron fluxes averaged over the 

plane of scan, vs time. Mean directional intensities 

are also indicated. The vertical arrows indicate the 

time of maximum for 158-keY and 79-keV electron 
♦ 

fluxes predicted by the theory of coherent propagation. 
The energies indicated are the appropriate mean energies 
for the channels . 

Fig* 7 Ji the directional proton fluxes averaged over the plane 
of scan, vs time. The dashed curves are constructed ■ 
using the curves in Figure 9 and indicate diffusive 

propagation. The onset of 1 -MeV protons is between 1600 

< ■ ) . ' 

and 2000 UT . The vertical arrows indicate the time of 
max i mum p re d i c t e d by the theory of coherent propagation. 
The energies indicated are the appropriate mean energies 
for the channels. 

Fig. 8 3 = v/c vs the time interval between the onset of the 

flare at the sun and the onset of the particle event in 
the different energy channels. 3 is calculated for the 
upper energy limit in each channel, i.e. P6 (46 MeV), 

P5 (13.3 MeV), P4 (5.4 MeV), P3 (1.35 MeV), El (102 keV) 
and E2 (185 keV). Channel 01 responds to protons with 
energies >100 MeV as indicated by a vertical line. The 
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• horizontal lines indicate the uncertainty in determining 

the onset times. 

Fig. 9 Test for diffusive propagation in the ADB-model. 0 is 

the average of the di recti onal , proton fluxes in the plane 
of scan. The vertical arrows indicate the points which 
correspond to the observed times of maximum. Only two 
arrows are shown because the diffusion maximum in the 
flux of 3-MeV protons, as evident from Figure 7, is not 
observed owing to the influence of the magnetic irregu- 
1 ari ty . 

Fig* 10 Directional proton fluxes vs cos 9* Here 0 is tfie angle 

mag nesM as a* a Vascos ano $ne 

tions in which the flux is measured. The distributions 
in front of the shock are shown. Satisfactory statisti- 
cal precision was 1 achieved 'by averaging the particle 
measurements over the time intervals indicated .in the 
figures . 

Fig. 11 A comparison with proton measurements from other sources 
during the event. The data points and curve 1 are from 
Durney et al. (1972), and represent the angular distri- 
bution in interpl anetary space of/vlOO-MeV protons at 
1115 UT on November 18, 1968, as determined from obser- 
vations of polar-cap proton fluxes together with tra- 
jectory calculations. Curves 2 and 3 are the distribu- 
tions of 25-MeV protons observed by 0G0-5 in interplane- ■ 
tary space at 1210 UT and 1340 UT, respectively* 
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F "f g . 12 Intensity vs cos 3 at a distance of ^300 in front 

of the propagating interplanetary shock. The data show 

a minimum in the angular distribution of 25-MeV protons 

in a direction nearly perpendicular to the direction of 

the mean magnetic field. 

Fig. 13 Proton angular distributions behind the shock - (a) well 
away from the magnetopause and (b) as the satellite 
approaches the magnetopause. 

Fig. 14 Angular distributions of electrons in front of the 
shock (a) shows that the directional flux intensity 
is decreasing as cos 0 ■+ -1 at early times following 
onset, while at later times, as illustrated in (b) 
and (c), the intensity Is independent of direction 
when cos 0 < 0. At all times the distribution for 
cos 0 >0 is linear in cos 8. ' • ; 

- ( „ ’ 1 , . 1 v ’ 

Fig. 15 Time variation of the electron anisotropy <5 = CJ max m n * n ) / 
^max + ^mi n ^ * 

Fig. 16 Angular distributions of electrons observed behind the 
i nterpl anetary shock, A bi-directional field-aligned 
anisotropy is clearly displayed only at 79 keV. 
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